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Abstract Measurements are presented of differential cross- 
sections for top quark pair production in pp collisions at 
= 7 TeV relative to the total inclusive top quark pair produc- 
tion cross-section. A data sample of 2.05 fb^' recorded by 
the ATLAS detector at the Large Hadron Collider is used. 
Relative differential cross-sections are derived as a function 
of the invariant mass, the transverse momentum and the ra- 
pidity of the top quark pair system. Events are selected in the 
lepton (electron or muon) + jets channel. The background- 
subtracted differential distributions are corrected for detec- 
tor effects, normalized to the total inclusive top quark pair 
production cross-section and compared to theoretical pre- 
dictions. The measurement uncertainties range typically be- 
tween 10% and 20% and are generally dominated by sys- 
tematic effects. No significant deviations from the Standard 
Model expectations are observed. 

1 Introduction 

The top quark is the most massive known fundamental con- 
stituent of matter Its unexplained large mass suggests an 
important connection to the electroweak symmetry breaking 
mechanism. The measurement of the top-antitop (ff) quark 
production cross-section (oi,) in various decay channels al- 
lows a precision test of perturbative QCD. In addition, the 
ff production process is an important background for Stan- 
dard Model (SM) Higgs boson searches, and in searches for 
physics beyond the SM. Also, a rich set of possible new 
particles and interactions might appear at the Large Hadron 
Collider (LHC) and modify the production and/or decay of 
top quarks. 

The inclusive f f production cross-section has been mea- 
sured by the ATLAS and CMS Collaborations with increas- 
mg precision in a variety of channels using data col- 
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lected in 20 1 and 20 1 1 . The unprecedented number of avail- 
able f f events (tens of thousands) enables detailed investiga- 
tions of the properties of top quark production in terms of 
characteristic variables of the f f system. This paper focuses 
on three observables of the ff system: the invariant mass 
(nitf), the transverse momentum (pT,rf) and the rapidity (ytf). 
To enable direct comparisons to theoretical models the dif- 
ferential distributions are unfolded for detector effects and 
corrected for acceptance effects. Theoretical predictions for 
the f f invariant mass distribution accurate to next-to-next-to- 
leading logarithm (NNLL) and next-to-leading order (NLO) 
are currently available |5], with a typical uncertainty of around 
12% at nift ~ 1 TeV. Comparisons of mass, transverse mo- 
mentum, and rapidity distributions are also made between 
unfolded data and NLO predictions taken from the MCFM 
generator In addition, the data are compared to predic- 
tions from the MC@NLO 101 and ALPGEN H generators 
with particular choices of parameter settings. 

The nitt distribution is sensitive to particles beyond the 
SM, such as new i-channel resonances that can modify the 
shape of the differential production cross-section in different 
ways depending on their spin and colour properties |9]. In 
addition to Tevatron experiment searches, both the ATLAS 
and CMS Collaborations have performed direct searches for 
specific narrow and wide resonances that extend mass lim- 
its to the TeV region |10-12]. The CDF Collaboration has 
performed a measurement of the differential cross-section 
as a function of m„- using the data collected in proton- 
antiproton (pp) collisions at ^/s = 1.96 TeV. The result is 
consistent with the SM expectation as predicted by PYTHIA 
(version 6.216) [ 14]. A potentially intriguing deviation from 
the SM prediction is observed in the measured forward-back- 
ward angular asymmetry between f and f quarks produced 
together in pp collisions at the Tevatron, particularly in events 
with large m„- [15]. Nearly all new physics scenarios that 
could explain this deviation should be observable at the LHC 
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as a resonant or non-resonant enhancement with respect to 
the SM in tt production at large m,f lll6ll . 



2 Detector, data and simulation samples 

The ATLAS detector iH at the LHC covers nearly the en- 
tire solid angle around the collision point. It consists of an 
inner tracking detector (ID) comprising a silicon pixel de- 
tector, a silicon microstrip detector, and a transition radia- 
tion tracker, providing tracking capability within pseudora- 
pidity Q 1 77 1 < 2.5. The ID is surrounded by a thin super- 
conducting solenoid providing a 2 T axial magnetic field, 
and by liquid argon (LAr) electromagnetic (EM) sampling 
calorimeters with high granularity. An iron/scintillator tile 
calorimeter provides hadronic energy measurements in the 
central pseudorapidity range (|t7| < 1.7). The end-cap and 
forward regions are instrumented with LAr calorimeters for 
both electromagnetic and hadronic energy measurements up 
to |t7| < 4.9. The calorimeter system is surrounded by a 
muon spectrometer incorporating three superconducting toroid 
magnet assemblies. 

A three-level trigger system is used to select high-/?T 
events. The level- 1 trigger is implemented in hardware and 
uses a subset of the detector information to reduce the rate 
to a design value of at most 75 kHz. This is followed by 
two software-based trigger levels, which together reduce the 
event rate to about 300 Hz. This analysis uses LHC pp col- 
lisions collected by the ATLAS detector between March and 
August 2011, corresponding to an integrated luminosity of 
2.05 fb-^ 

Simulated top quark pair events are generated using the 
MC@NLO Monte Carlo (MC) generator version 3.41 with 
the NLO parton distribution function (PDF) set CTEQ6.6 1 18], 
where the top mass is set to 172.5 GeV. Renormalization and 
factorization scales are set to the same value: the square root 
of the average of the t and t quarks squared transverse en- 
ergies. Parton showering and the underlying event are mod- 
eUed using HERWIG H and JIMMY Hi using the AUETl 
tune 12^1], respectively. The tt sample is normalized to a 
cross-section of 164.6 pb, obtained with an approximate NNLO 
prediction 11221] . Single top events are also generated using 
MC@NLO, while the production of W /Z bosons in asso- 
ciation with jets is simulated using the ALPGEN genera- 
tor interfaced to HERWIG and JIMMY with CTEQ6L1 



^ ATLAS uses a right-handed coordinate system with its origin at the 
nominal interaction point in the centre of the detector and the z-axis 
along the beam pipe. The x-axis points to the centre of the LHC ring, 
and the y-axis points upward. Cylindrical coordinates {r,(j>) are used 
in the transverse plane, being the azimuthal angle around the beam 
pipe. The pseudorapidity is defined in terms of the polar angle 6 as 
77 = — lntan(6/2). Transverse momentum and energy are defined as 
Pt = p sin 9 and Ej = EsinG, respectively. The distance AR is defined 
as AR = ^/ {A(j))^ + {ArjY, where A(j) and Ar} are the separation in 
azimuthal angle and pseudorapidity, respectively. 



PDFs ll23h . W+jets events containing bb pairs, cc pairs and 
single c-quark (heavy flavour) were generated separately us- 
ing matrix elements with massive b- and c-quarks. An overlap- 
removal procedure is used to avoid double counting due to 
heavy quarks from the parton shower Diboson events (WW, 
WZ, ZZ)ave generated using HERWIG with MRST LO* 

PDFs ilii . 

All Monte Carlo simulation samples are generated with 
multiple pp interactions per bunch crossing (pile-up). These 
simulated events are re-weighted so that the distribution of 
the average number of interactions per pp bunch crossing 
in simulation matches that observed in the data. This aver- 
age number varies between data-taking periods and ranges 
typically between 4 and 8. The samples are then processed 
through the GEANT4 f23\ simulation of the ATLAS detec- 
tor \2M and the standard ATLAS reconstruction software. 



3 Event selection 

Events are selected in the lepton (electron or muon) + jets 
channel. The reconstruction of tt events in the detector is 
based on the identification and reconstruction of electrons, 
muons, jets and missing transverse momentum. The defini- 
tions of these objects are identical to those used in Ref. 1271] . 
The same event selection as in Ref. [27J is used with the ad- 
dition of a requirement on the kinematic likelihood resulting 
from the event reconstruction described in Sect.|5j 



3.1 Object definitions 

Electron candidates are defined as energy deposits in the 
EM calorimeter associated with well-reconstructed tracks of 
charged particles in the ID. The candidates are required to 
meet stringent identification criteria based on EM shower 
shape information, track quality variables and information 
from the transition radiation tracker |28]. All candidates are 
required to have Ej > 25 GeV and |rjciu| < 2.47, where Tjciu 
is the pseudorapidity of the EM calorimeter cluster associ- 
ated with the electron. Candidates in the transition region 
between the barrel and end-cap calorimeters 1.37 < |t7c1u| < 
1.52 are rejected. 

Muon candidates are reconstructed by combining track 
segments in different layers of the muon chambers. Such 
segments are assembled starting from the outermost layer, 
with a procedure that takes material effects into account, and 
are then matched with tracks found in the ID. The candidates 
are then re-fitted, exploiting the full track information from 
both the muon spectrometer and the ID, and are required to 
have pt > 20 GeV and |tj | < 2.5. 

Jets are reconstructed with the anti-A:, algorithm 1291] 



with a distance parameter of 0.4 using clusters formed from 
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calorimeter cells with significant energy deposits ("topoclus- 
ters") at the EM scale. The jet energy is then corrected to the 
hadronic scale using pj- and Tj -dependent correction factors 
derived from simulation and validated with data llsoll . 

The missing transverse momentum and its magnitude 
is™^** are derived from topoclusters at the EM scale and cor- 
rected on the basis of the energy scale of the associated 
physics object, if any |31]. Contributions from muons are 
included using their momentum measured from the tracking 
and muon spectrometer systems. The remaining clusters not 
associated with high-/7T objects are added at the EM scale. 

Both the electron and muon candidates are required to 
be isolated to reduce the backgrounds from hadrons mim- 
icking lepton signatures and leptons from heavy-flavour de- 
cays. For electron candidates, the total transverse energy de- 
posited in the calorimeter in a cone of AR = 0.2 around the 
electron candidate is required to be less than 3.5 GeV after 
correcting for the energy associated with the electron and 
for energy deposited by pile-up. For muon candidates, the 
isolation is defined in a cone of AR = 0.3 around the muon 
direction. In that region both the sum of track transverse mo- 
menta for tracks with pj > 1 GeV and the total energy de- 
posited in the calorimeter are required to be less than 4 GeV, 
after subtracting the contributions from the muon itself. 

Jets within AR = 0.2 of an electron candidate are re- 
moved to avoid double counting electrons as jets. Subse- 
quently, muons within AR = 0.4 of the center of a jet with 
Pj > 20 GeV are removed in order to reduce the contamina- 
tion caused by muons from hadron decays. 

The reconstruction of tt events is aided by the abiUty to 
tag jets from the hadronization of ^-quarks using the com- 
bination of two /j-tagging algorithms 1321]. One fe-tagger de- 
rives the properties of vertices related to b- and c-hadron 
decays inside jets by assuming the vertices to lie on a line 
connecting them to the primary vertex H. A likelihood dis- 
criminant between b-, c- and light-quark jets is derived by 
using the number, the masses, the track energy fraction, the 
flight-length significances and the track multiplicities of the 
reconstructed vertices as inputs. The other ^-tagging algo- 
rithm employs the transverse and longitudinal impact pa- 
rameter significances of each track within the jet to derive 
a likelihood that the jet originates from a b-qnaik. The re- 
sults of the two taggers are combined, using a neural net- 
work, into a single discriminating variable. The combined 
tagger operating point chosen for the present analysis cor- 
responds to a 70% tagging efficiency for b-jets in simulated 
tt events, while light-flavour jets (c-jets) are suppressed by 
approximately a factor of 100 (5). 

primary vertex is defined as a vertex reconstructed from a number 
of high-quality tracks such that the vertex is spatially compatible with 
the luminous region of interaction. Primary vertices in an event are 
ordered by Xtrk Pj n-i;, where px.trk is the transverse momentum of an 
associated track. 



3.2 Selection of tt candidates 

The lepton + jets channel selection requires the appropri- 
ate single-electron or single-muon trigger to have fired (with 
thresholds at 20 GeV and 18 GeV respectively). Events pass- 
ing the trigger selection are required to contain exactly one 
reconstructed electron (muon) with Ej > 25 GeV {pj > 20 
GeV). The events are required to have at least one recon- 
structed primary vertex. The primary vertex, corresponding 
to that with highest Z^trk Pj i± is required to be reconstructed 
from at least five high-quality tracks. Jet quality criteria are 
applied to the data and events are discarded if any jet with 
Pj > 20 GeV is identified to be due to calorimeter noise or 
activity out of time with respect to the LHC beam cross- 
ings l3Qi1 . The E™^^ is required to be larger than 20 (35) 
GeV in the jj. -l-jets (e + jets) channel. The W boson trans- 
verse mass (nij), derived from the lepton transverse mo- 
mentum and the £'™'*'* |33], is required to be larger than 60 
GeV - fi^'"" (25 GeV) in the fi + jets (e + jets) channel. The 
requirements for the e+jets channel is more stringent in 
order to reduce the larger fake-lepton background. Events 
are required to have at least four jets with pj > 25 GeV 
and \ri\ < 2.5, where at least one of these jets is required 
to be ^-tagged. Finally, events are retained only if they have 
a kinematic likelihood ln(L) > —52 resulting from the event 
reconstruction described in Sect.|5]. 



4 Background determination 

The main expected backgrounds in the lepton + jets channel 
are W+jets which can give rise to the same final state as the 
tt signal, and fake leptons. They are both estimated using 
auxiliary measurements. The other backgrounds are of elec- 
troweak origin and are estimated from simulation. All back- 
ground determination methods are identical to those used in 
Ref. Q. 



4. 1 Fake-lepton background 

The multijet background with misidentified and non-prompt 
leptons (referred to collectively as fake leptons) in both the 
e + jets and ji + jets channels is evaluated with a matrix 
method, which relies on defining loose and tight lepton sam- 
ples m [Istl and measuring the fractions of real (Ereai) and 
fake (Efake) loose leptons that are selected as tight leptons. 
The fraction Ereai is measured using data control samples 
of Z boson decays to two leptons, while gfake is measured 
from data control regions dominated by the contributions of 
fake leptons. Contributions from W+jets and Z+jets back- 
grounds are subtracted in the control regions using Monte 
Carlo simulation. 
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For the jj, +jets channel, the loose data sample is defined 
by discarding the isolation requirements in the default muon 
selection. The fake-muon efficiencies are derived from a low- 
control region, irij < 20 GeV, with an additional re- 
quirement /i™^^ + < 60 GeV. The efficiencies for real 
and fake muons are parameterized as a function of muon | rj \ 
and of the leading jet p-^. 

For the e + jets channel, the loose data sample is defined 
by selecting events with electrons meeting looser identifica- 
tion criteria. The 3.5 GeV electron isolation requirement is 
loosened to 6 GeV. The fake-electron efficiencies are deter- 
mined using a low-is™^^ control region (5 GeV < E™*^** < 
20 GeV). The efficiencies for real and fake-electrons are pa- 
rameterized as a function of electron | T] | . 



4.2 W+jets background estimation 

The W+jets background estimation consists of three steps. 

The first step is to determine the flavour composition 
of the W+jets background in the signal region before b- 
tagging. Since the theoretical prediction for heavy flavour 
fractions in W+jets suffers from large uncertainties, a data- 
driven approach was developed to constrain these fractions 
with inputs from MC simulation. Samples with a lower jet 
multiplicity, obtained from the selection described in Sect. 13. 21 
but requiring exactly one or two jets instead of four or more 
jets, are analysed. 

The numbers Wj.^p^tta^' ^/'tagged' of ^ + '"jets events in 
these samples (with i—l, 2), before and after applying the b- 
tagging requirement, are calculated from the observed events 
by subtracting the small contributions from other Standard 
Model processes — electroweak (WW, WZ, ZZ, and Z+jets) 
and top (ff and single top) processes — predicted by the 
simulation and by subtracting the fake-lepton background 
obtained as described in Sect. 14.11 

A system of three equations — expressing the number of 
W+l-jet events after /^-tagging and W+2-jets events before 
and after Zj-tagging — can be written with eight independent 
flavour fractions as the unknowns, corresponding to frac- 
tions of Wi>^+jets, Wcc +jets, Wc +jets and W+light-jets 
events in the one- and two- jet bins before /^-tagging. In the 
equations involving tagged events, the simulation prediction 
is used to include the eight tagging probabilities of the dif- 
ferent W+jets event types. For each flavour, the fractions 
in the one-jet and two-jet bins are related using the simula- 
tion's prediction of their ratio. These predictions reduce the 
number of independent fractions to four Finally, the ratio of 
the Wcc +jets to the Wbh+jets fractions in the two-jet bin is 
fixed to the value obtained from simulated events in order to 
obtain three independent fractions in the three equations. 

The resulting scale factors for the heavy flavour fractions 
in simulated W+jets events are 1.63 ±0.76 for Wbh+jets 



and Wcc+jets events and 1.11 ±0.35 for Wc+jets events. 
These are applied to the relevant Monte Carlo samples. The 
uncertanties on these scale factors are derived from system- 
atic variations of the inputs to the method (see Sect. I6.2l i. 
The fraction of W + light-jets events is scaled by a factor 
0.83 to keep the total number of pre-tagged Monte Carlo 
W+jets events fixed. When applied to the signal region, an 
additional 25% uncertainty is applied to these fractions, cor- 
responding to the uncertainty in the Monte Carlo prediction 
for the ratio of flavour fractions in different jet multiplicities. 

The second step is to determine the overall normaliza- 
tion of W+jets background in events with four or more jets 
before ^j-tagging. At the LHC the rate of W++jets events is 
larger than that of W^ +jets events because there are more 
up-type valence quarks in the proton than down-type valence 
quarks. The ratio of W++jets to W^+jets cross-sections is 
predicted much more precisely than the total W+jets cross- 
section fi^-'s^. This asymmetry is used to measure the to- 
tal W+jets background from the data. To a good approxi- 
mation, processes other than W+jets give equal numbers of 
positively and negatively charged leptons. Consequently the 
total number of W+jets events in the selected sample can be 
estimated as 



W>4,pre-tag = M 



w- 



f ^MC + 1 



D 



(1) 



The charge-asymmetric single top contribution is estimated 
from simulation and subtracted. The values D+ {D^ ) are the 
total numbers of events in data meeting the selection criteria 
described in Sect. 13.21 before the /j-tagging and likelihood 
requirement, with positively (negatively) charged leptons. 

The value of rivir = ^S^^^^w-Tvl is derived from Monte 
Carlo simulation, using the same event selection. The ratio 
rMC is 1.56 + 0.06 in the e+jets channel and 1.65 + 0.08 
in the /i +jets channel. The largest uncertainties on r^c de- 
rive from uncertainties in PDFs, the jet energy scale, and the 
heavy flavour fractions in W+jets events. 

Finally, in the third step, the number of W +jets events 
passing the selection with at least one i>-tagged jet is deter- 
mined to be 13311 



^>4,tagged — W^>4,pre-tag ' /2,tagged ' ^2^>4- 



(2) 



The value ^.tagged = W^%^JW^^,,^„ is the fraction of 



— lydata 

2, tagged/ ''''2, pre-tag 

W+ 2 jets events meeting the requirement of having at least 
one Z7-tagged jet, and ^2^>4 = />4^tagged//2^tagged is the ratio 
of the fractions of simulated W +jets events passing the re- 
quirement of at least one /7-tagged jet, for at least four and 
exactly two jets, respectively. The value of /2.tagged is found 
to be 0.063 ± 0.005 in the e + jets channel and 0.068 ± 0.005 
in the pL +jets channel. The ratio k2^>4 is found to be 2.52 ± 
0.36 in the e + jets channel and 2.35 ± 0.34 in the ji + jets 
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channel. The uncertainties include both systematic contri- 
butions and contributions arising from the limited number 
of simulated events. 



4.3 Other backgrounds 

The numbers of background events from single top produc- 
tion, Z+jets and diboson events are evaluated using Monte 
Carlo simulation. The prediction for Z+jets events are nor- 
malized to the approximate NNLO cross-sections as deter- 
mined by the FEWZ program 1 37 ] , using the MSTW2008NLO 
PDFs 1,34. 38]. The prediction for diboson events is normal- 
ized to the NLO cross-section as determined by the MCFM 
program fs?] using the MSTW2008NLO PDFs. The approx- 
imate NNLO cross-section results from Refs. 1 43-42] are 
used to normalize the predictions for single top events. 



5 Reconstruction 

Measurements of differential cross-sections in top quark pair 
events require full kinematic reconstruction of the tt system. 
The reconstruction is performed using a likelihood fit of the 
measured objects to a theoretical leading-order representa- 
tion of the tt decay. The same reconstruction method as in 
Ref. II27II is used. The likelihood is the product of three fac- 
tors. The first factor is the product of Breit-Wigner distribu- 
tions for the production of W bosons and top quarks, given 
the four-momenta of the true tt decay products. The sec- 
ond factor is the product of transfer functions representing 
the probabilities for the given true energies of the tt decay 
products to be observed as the energies of reconstructed jets, 
leptons and as missing transverse energy. The third factor is 
the probability to b-tag a certain jet, given the parton it is 
associated with. The pole masses of the W bosons and the 
top quarks in the Breit-Wigner distributions are set to 80.4 
GeV and 172.5 GeV, respectively. 

The likelihood is maximized by varying the energies of 
the partons, the energy of the charged lepton, and the com- 
ponents of the neutrino three-momentum. The maximization 
is performed over all possible assignments of jets to partons, 
and the assignment with the largest likelihood is used for all 
further studies. The four-momenta of the top quarks are then 
obtained by summing the four momenta of the decay prod- 
ucts resulting from the kinematic fit. The unconstrained z 
component of the neutrino momentum is a free parameter in 
the fit. 

Simulation studies aimed at enhancing the fraction of 
properly reconstructed tt events are used to determine a re- 
quirement on the likelihood of the kinematic fit. The like- 
lihood distribution for the events after selection, except for 
the likelihood requirement ln(L) > —52, is shown in Fig.[T] 



The likelihood optimally encapsulates all relevant informa- 
tion about the data agreement with simulation. Distributions 
of the reconstructed invariant mass, transverse momentum 
and rapidity of the reconstructed top-antitop pair, after all 
selection requirements, are shown in Fig.|2] 

The numbers of expected and observed data events in 
each channels after pre-tag, tagged and full event selection 
are listed in Table [1] The data agrees with the expectation 
within the systematic uncertainties. 



6 Systematic uncertainties 

For each systematic effect the analysis is re-run with the 
variation corresponding to the one standard deviation change 
in each bin. The varied distributions are obtained for the up- 
ward and downward shift for each effect, and for each chan- 
nel separately. If the direction of the variation is not defined 
(as in the case of the estimate resulting from the difference 
of two models), the estimated variation is assumed to be the 
same size in the upward and the downward direction and is 
symmetrized. The baseline distribution and the shifted dis- 
tributions are the input to the pseudo-experiment calculation 
(see Sect. |8]l that performs unfolding, efficiency correction, 
and enables combination of the e +jets and ji + jets chan- 
nels. 

The sources of systematic uncertainties are arranged in 
approximately independent groups. They are further catego- 
rized into detector modelling, and modelling of signal and 
background processes. The estimation of the variations re- 
sulting from the systematic uncertainty sources is the same 
as Ref. 127J. 



6.1 Detector modelling 

Muon and electron trigger, reconstruction and selection effi- 
ciencies are measured in data using Z and W decays and in- 
corporated into the simulation using weighted events. Each 
simulated event is weighted with the appropriate ratio (scale 
factor) of the measured efficiency to the simulated one. The 
uncertainties on the scale factors are estimated by varying 
the lepton and signal selections and background uncertain- 
ties. For lepton triggers the systematic uncertainties are about 
1 %. The same procedure is used for lepton momentum scale 
and resolution scale factors resulting in uncertainties of 1- 
1.5%. The corresponding scale factor uncertainties for elec- 
tron (muon) reconstruction and identification efficiency are 
1% (0.5%) and 3% (0.5%) respectively. 

Information collected from collision data, test-beam data, 
and simulation is used to determine the jet energy scale; its 
uncertainty ranges from2.5% to 8%, varying with jet pj and 
rj 1 30]. The uncertainties include flavour composition of the 
sample and mis-measurements due to nearby jets. Pile-up 
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Fig. 1 Negative logarithm of the Hkehhood obtained from the kinematic fit described in the text. Data are compared to expectation from Monte 
Carlo simulation and data-driven expectation for the (a) e + jets and (b) fi +jets selected events. All selection criteria are applied except for the 
likelihood requirement. The band represents the 68% confidence level inteiTal of total uncertainty on the prediction. The bin corresponding to the 
largest — In(likelihood) value includes events with — In(likelihood) > 70 and the associated prediction. 



Table 1 Numbers of predicted and observed events. The selection is shown after applying pre-tag, tagged, and the full selection criteria including 
the likelihood requirement. The quoted uncertainties include all uncertainties given in Sect. |6] except those from PDF and theory normalization. 
The numbers correspond to an integrated luminosity of 2.05 fb^' in both e+jets and ^ +jets samples. 



Channel 


/X + jets pre-tag 


/i + jets tagged 


/i + jets L-req 


e + jets pre-tag 


e+jets tagged 


e+jets L-req 




tt 


15800 ± 1300 


13900+ 1100 


11100 + 700 


10700 + 900 


9400 + 800 


7400 + 500 


W+jets 


19000 ± 5000 


3000 ± 1200 


1700 + 700 


13000 ± 3300 


2200 + 900 


1300 + 500 


Single top 


950 ± 70 


760 ± 80 


490 + 50 


660 ± 50 


530 + 50 


338 + 32 


Z+jets 


2200 ± 200 


309 ± 34 


192 + 20 


1750 + 330 


240 + 50 


154 + 26 


Diboson 


298 ± 28 


53 + 7 


34 + 4 


181 + 19 


32 + 5 


21 + 3 


Fake-leptons 


3400 ± 1700 


1100+ 1100 


800 + 800 


2000 + 1000 


400 + 400 


250 + 250 


Signal+bkg 


42000 ± 6000 


19200 + 2600 


14400 + 1700 


28000 + 4000 


12800 + 1700 


9500+ 1100 


Observed 


42327 


19254 


14416 


26488 


12457 


9187 



gives an additional uncertainty of 2.5% (5%) in the cen- 
tral (forward) region. An extra uncertainty of up to 2.5% 
is added to account for the fragmentation of ^j-quarks. The 
jet energy resolution and reconstruction efficiency are mea- 
sured in data using the same methods as in Refs. I sol 431. 



Jet energy resolution uncertainties range from 9-17% for jet 
Pt ~ 30 GeV to about 5-9% for jet pj > 180 GeV depend- 
ing on jet J]. The jet reconstruction efficiency uncertainty is 
1-2%. The uncertainties from the energy scale and resolu- 
tion corrections on Ieptons and jets are propagated to the un- 
certainties on missing transverse momentum. Uncertainties 
on Ej'^^ also include contributions arising from calorime- 
ter cells not associated to jets and from soft jets (those in the 
range 7 GeV < pj <20 GeV). The /^-tagging efficiency scale 



factors have uncertainties between 6% to 15%, and mis-tag 
rate scale factor uncertainties range from 10% to 21%. The 
scale factors are derived from data and parameterized as a 
function of jet pj. 

A small region of the liquid argon calorimeter could not 
be read out in a subset of the data corresponding to 42% of 
the total dataset. Corresponding data and simulated events 
where a jet with pj > 20 GeV is close to the failing re- 
gion are rejected. This requirement rejects about 6% of the 
events. A systematic uncertainty is derived from variations 
of the /?x-threshold of the jets by 20% resulting from studies 
of the response of jets close to the failing region, using dijet 
p-Y balance in data. 
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Fig. 2 Distributions of the reconstructed (a-b) tt mass, rriff, (c-d) the transverse momentum, px,?/. (e-f) the //"rapidity, >Vf, before background 
subtraction and unfolding. Data are compared to the expectation derived from simulation and data-driven estimates. All selection criteria are 
applied for the (a, c, e) e + jets and (b, d, f) /i H-jets channels. The uncertainty bands include all contributions given in Sect. [6l except those from 
PDF and theory normalization. 
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The uncertainty on the measured luminosity is 3.7% 144 



6.2 Signal and background modelling 

Sources of systematic uncertainty for the signal are the choice 
of generator, parton shower model, hadronization and un- 
derlying event model, the choice of PDF, and the tuning of 
initial- and final-state radiation. Predictions from the MC- 
@NLO and POWHEG Ii46|,i47j generators ai-e compared 
to determine the generator dependence. The parton shower- 
ing is assessed by comparing POWHEG samples interfaced 
to HERWIG and PYTHIA, respectively. The amount of 
initial- and final-state radiation is varied by modifying pa- 
rameters in AcerMC [48] interfaced to PYTHIA. The pa- 
rameters are varied in a range comparable to those used in 
the Perugia Soft/Hard tune variations [49]. The impact of the 
PDF uncertainties is studied using the procedure described 
inRef. (sd]. 

Background processes are either estimated by simula- 
tion or using auxiliary measurements, see Sect.|4] The uncer- 
tainty on the fake-lepton background is estimated by vary- 
ing the requirements on the low-m^ and low-£'™''*'* control 
regions, taking into account the statistical uncertainty and 
background corrections. The total uncertainty is estimated 
to be 100%. The normalization of W+jets is derived from 
auxiliary measurements using the asymmetric production of 
positively and negatively charged W bosons in W+jets events. 
The uncertainties are estimated to be 21% and 23% in the 
four-jet bin, for the electron and muon channels respectively. 
This uncertainty is estimated by evaluating the effect on both 
rue and fc2^>4 from the jet energy scale uncertainty and dif- 
ferent PDF and generator choices. Systematic uncertainties 
on the shape of W^-|-jets distributions are assigned based on 
differences in simulated events generated with different fac- 
torization and parton matching scales. Scaling factors cor- 
recting the fraction of heavy flavour contributions in simu- 
lated W+jets samples are derived from auxiliary measure- 
ments (see Sect. l4.2l i. The systematic uncertainties are found 
by changing the normalizations of the non-W processes within 
their uncertainties when computing W,-^p'g_(2g and Wf^^^^^^, 
as well as taking into account the impact of uncertainties 
in i>-tagging efficiencies. The total uncertainties are 47% for 
Wbh+jets and Wcc+jets contributions and 32% for Wc+jets 
contributions. In the ji +jets channel the fractional contri- 
butions of Wbb+jets, Wcc +jets and Wc+jets samples to 
the total W+jets prediction are estimated to be 9%, 17% 
and 12% (36%, 25% and 17%) respectively, before (after) 
the ^-tagging requirement. In the e + jets channel the frac- 
tional contributions of Wbb+jets, Wcc +jets and Wc+jets 
samples to the total W+jets prediction are estimated to be 
9%, 17% and 13% (35%, 25% and 17%) respectively, be- 
fore (after) the i>-tagging requirement. The normalization of 



Z+jets events is estimated using Berends-Giele-scaling 15 111 . 
The uncertainty in the normalization is 48% in the four-jet 
bin and increases with the jet multiplicity. The uncertainties 
on the normalization of the small background contributions 
from diboson and single top production are estimated to be 
about 5% [34,18, 39J and 10% L4fl-i42J, respectively. 

The statistical uncertainty on the Monte Carlo predic- 
tion due to limited Monte Carlo sample size is included as a 
systematic uncertainty in each bin for each process. 



7 Cross-section unfolding 

7.1 Unfolding procedure 

The underlying binned true differential cross-section distri- 
butions (Oj) are obtained from the reconstructed events us- 
ing an unfolding technique that corrects for detector effects. 
The unfolding starts from the reconstructed event distribu- 
tion (A^,), where the backgrounds (B,) have been subtracted. 
The unfolding uses a response matrix (^,7), see Eq. (O, de- 
rived from simulated tt events, which maps the binned gen- 
erated events to the binned reconstructed events. The kine- 
matic properties of the generated t and f partons in simulated 
tt events define the "true" properties of the tT events. In its 
simplest form the unfolding equation can be written as 



j j 



(3) 



where ^ is the integrated luminosity, Mjj is the bin migra- 
tion matrix (see Fig. O, and Aj is the acceptance which in- 
cludes the branching fraction BR{tt — >> 1/) = 0.543, taken 
from Ref. fs?], for tt events in which at least one W boson 
decays to leptons. 

The estimated acceptances for simulated tT events as a 
function of m„-, pj fj and y,,- are reported in Table |2] The 
overall acceptances before the requirement on the likelihood 
value are comparable to previous measurements llisll . The 
additional requirement on the likelihood value is expected 
to retain a large fraction of the previously selected tt events 
(see Table[T]i. A finely binned illustration of the acceptances 
is shown in Fig. ID The reduction in acceptance associated 
with high niif and pjjf values is predominantly due to the 
presence of increasingly non-isolated leptons coupled to lower 
jet multiplicity as tt decay products are forced in a closer 
space region by the boost at large top quark pj. In the case 
of high |y,f| it is mainly due to jets falling outside of the 
required pseudorapidity range (see Sect. 13.21 ). 

The cross-section (7, is then extracted by solving Eq. (O 



(4) 



The bin size is optimized using pseudo-experiments drawn 
from simulated events including systematic uncertainties. 
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Fig. 3 Migration matrices for (a-b) m,j, (c-d) pT,rf, and (e-f) y,j estimated from simulated tt events passing all (left)e+jets and (right) ju -I- jets 
selection criteria. The unit of the matrix elements is the probabiUty for an event generated at a given value to be reconstructed at another value. 
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Table 2 Table of acceptances for m,f, pj ,, and y,,-. The acceptance is defined according to Eq.|3]and includes the branching fraction BR{tt — >> 
1/) = 0.543 from Ref. 1 52] for ti events in which at least one W boson decays to leptons. In the case of acceptances in positive and negative 
symmetric bins are consistent within uncertainties. 



m„- [GeV] Acceptance [%] 
e+jets /i+jets 
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Fig. 4 Acceptance as a function of (a) tt mass, m,,-, (b) the tt transverse momentum, pjji, and (c) the tt rapidity, y,j. The acceptance is defined 
according to Eq.[3]and includes the branching fraction BR{tt — >> 1/) = 0.543 from Ref. |52] for tt events in which at least one W boson decays 
to leptons. The error bars show only the uncertainty due to limited Monte Carlo sample size. 
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The adopted optimization strategy is to choose as small a 
bin size as possible without substantially increasing the total 
uncertainty after unfolding. This effectively means keeping 
about 68% of the events on the diagonal of the migration ma- 
trix, and requiring that the condition numbe:[| of the migra- 
tion matrix is 0{l). The finely binned distributions before 
unfolding reported in Fig. |2] show good agreement between 
reconstructed data and the MC and data-driven predictions. 

To evaluate the performance of the unfolding procedure, 
and to estimate the systematic uncertainties, Eq. (|4]i has been 
extended to the following form to allow detailed studies us- 
ing pseudo-experiments 



(7j(4) 



-^)jiidk)[P{N,)-B,idk)] 



Aj{dk)^{dk) 



(5) 



where P{Ni) is the Poisson distribution with mean A^,, and 
dk are continuous variables representing the systematic un- 
certainties, drawn from a Gaussian distribution with zero 
mean and unit standard deviation. A cross-section estimate 
Oj is extracted for a given variable {m,t, PT,tt, yrf) from each 
pseudo-experiment. The distribution of (7, resulting from 
the pseudo-experiments is an estimator of the probability 
density of all possible outcomes of the measurement. Two 
thousand pseudo-experiments are used to extract the cross- 
section values. The 68% confidence interval provides the 
cross-section uncertainty. The parametric dependence on d^ 
in {M^^)ij, and other functions, is approximated using the 
linear term in the Taylor expansion, treating positive and 
negative derivative estimates separately. 

A closure test is performed by unfolding simulated (folded) 
events where c/^ = 0. The deviation of the unfolded cross- 
section from the known true cross-section input, used for the 
detector simulation folding, is consistent with zero within 
1% uncertainty. The most important test of the unfolding is 
to test the ability to unfold a distribution significantly dif- 
ferent from the Monte Carlo expectation. This is done by re- 
weighting simulated tt events so that the number of events in 
a single bin of true m,,- is doubled. The observed linearity of 
the response to these "delta-like" pulses is within 1%. The 
same test was also performed using a regularized unfolding 
technique based on Singular Value Decomposition jsH. The 
size of the "delta-like" pulses was then found to be substan- 
tially reduced (at least by 30%) after unfolding, even under 
the mildest regularization conditions. Given the bias from 
regularized unfolding, all final results are derived using the 
plain matrix inversion described above. The increased sta- 
tistical uncertainty of this unregularized result is tolerated 
given that the total uncertainty is dominated by systematic 
effects. 



^The condition number k is defined as ^ = | |M| | • | ' 1 1 , and is a mea- 
sure of how much the matrix inversion increases the size of the uncer- 
tainties in the error propagation. 



7.2 Combination of channels 

The unfolded cross-sections from the two channels, e + jets 
and jj. +jets, are combined using a weighted mean which 
includes the full covariance matrix between the channels. 
Since the covariance matrix is used in the weighting, the 
estimate is a best linear unbiased estimator of the cross- 
section. The covariance matrix is determined in simulated 
events using the same pseudo-experiment procedure outlined 
in the previous section and derived from Eq. (|5]l. 



8 Results 

To reduce systematic uncertainties only relative cross-sections 
(differential cross-section normalized to the measured inclu- 
sive cross-section) are reported. The full procedure for the 
differential measurement is also contracted down to one bin 
to perform the measurement of the inclusive cross-section 
by using Eq. Q and Eq. (|4|. In this case the measurement is 
reduced to a standard "cut-and-count" technique (as used for 
the first ATLAS tt cross-section measurement |33]) and the 
response matrix is reduced to the standard acceptance cor- 
rection. The total inclusive cross-section, combining e + jets 
and ji +jets channels, is found to be (7„- = 160 ± 25 pb. 
The quoted uncertainty includes both statistical and system- 
atic contributions and the result is compatible with the ex- 
pected ff inclusive cross-section and with previous measure- 
ments 

The relative differential cross-section results are listed 
in Table [3] as a function of m,,-, pj ,f and yff. Both single- 
channel results and results from the combination are shown. 
A graphical representation for the combined results is shown 
in Fig.|5] The measurements are reported with their full un- 
certainty, combining statistical and systematic effects, and 
they are compared to NLO predictions from MCFM [6] for 
all variables; NLO-I-NNLL predictions from Ref. jstl are in- 
cluded for 1 /a da/dintf. Theory uncertainty bands include 
uncertainties on parton distribution functions, the strong cou- 
pling constant as and on factorization and renormalization 
scales. For the NLO predictions, the uncertainty on PDFs is 
set to the maximal spread of the predictions from three dif- 
ferent NLO PDF sets (CTEQ66, MSTW2008 an NNPDF) 
according to the PDF-specific recipe in Ref. ifsoll . Renormal- 
ization and factorization scales are set to the top quark mass 
value of 172.5 GeV and associated uncertainties are derived 
from an upward and downward scale variation of a factor of 
two. The overall NLO uncertainty is obtained by summing 
the contributions from PDF and scales in quadrature for vari- 
ations in the same direction. For the NLO-I-NNLL estimates 
the uncertainties are derived according to the approach of 
Ref. ij. The uncertainty on the MSTW2008NNLO PDFs 
and the strong coupling constant as at the 68% confidence 
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level is combined in quadrature with the uncertainties de- 
rived from the variations of the factorization scale and the 
renormalization scales. For 1/(7 do/dnitf the scale uncer- 
tainties are dominant. Predictions from MC @NLO and ALP- 
GEN are shown for fixed settings of the generators' parame- 
ters. The settings for MC @NLO are given in Sect.|51 ALP- 
GEN is version 2.13 using the CTEQ6L1 PDF with the top 
mass set to 172.5 GeV. Renormalization and factorization 
scales are set to the same value: the square root of the sum 
of the squared transverse energies of the final state partons. 
The matching parameters 1I54I1 for up to five extra partons are 
set to Ej"'^ = 20 GeV and Rmatch = 0.7. Parton showering and 
underlying event are simulated by HERWIG and JIMMY 
respectively, using the generator tune AUETl 112 ill . 

The impact of the different uncertainty sources on the 
final results is estimated and shown in Table |5] For 1 /a 
dc /diritt the relative statistical uncertainty varies from about 
2% at low m,f to about 20% at the highest m,,-, while the sys- 
tematic uncertainty ranges between 10% at intermediate m„- 
values to about 37% at the highest m,t. In relation to l/c7 
da /dp J , f the relative statistical uncertainty ranges between 
about 4% at low pjjf and about 12% at the highest pj ti val- 
ues, while the systematic uncertainty increases from about 
13% to 20% in the same interval. In the case of 1 /a da/dytf 
the relative statistical uncertainty increases from about 3% 
at low y,f to about 5% at the highest y,f values, while the 
systematic uncertainty changes from 4% to 10% over the 
same interval. Jet-related uncertainties are dominant for ot„- 
and PT/f, while for y,,- the dominant contributions are from 
fake-leptons and final-state radiation in addition to the jet 
uncertainties. 

No significant deviations from the SM expectations pro- 
vided by the different MC generators are observed. The SM 
prediction for the relative cross-section distribution can be 
tested against the measured values by using the covariance 
matrix between the measured bins. The correlation coeffi- 
cients are estimated using five thousand pseudo-experiments, 
see Table|4] The covariance matrices are derived by combin- 
ing the correlation coefficients with the uncertainties for the 
respective measurements reported in Table [3] for p^jf 
and ytf respectively. 

9 Conclusions 

Using a dataset of 2.05 fb^', the relative differential cross- 
section for tt production is measured as a function of three 
properties of the tt system: mass (m,,), pj ipj,tf) and rapid- 
ity (ytt)- The background-subtracted, detector-unfolded val- 
ues of 1/ff da/dniff, 1/(7 da/dpjtf and 1/(7 da/dy,, are 
reported together with their respective covariance matrices, 
and compared to theoretical calculations. The measurement 
uncertainties range typically between 10% and 20% and are 
generally dominated by systematic effects. No significant 



deviations from the SM expectations provided by the dif- 
ferent MC generators are observed. 
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ALPGEN are shown for fixed settings of the generators' parameters (details are found in Sect.[8j. 
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Table 3 Relative differential cross-section (top) l/ff d<j/dm,j, (middle) l/c da/dpj,f and (bottom) l/a da/dy,j measured in the e-|-jets, 
/I -hjets and the combined ^'-|-jets channel. 





[GeV] 




1/(7 da/dm,j [1/TeV] 








e-hjets 


^ +jets 


^+jets 




250 


-450 


2.2 ± 0.4 


2.5 -fO.3 / -0.4 


2.4 +03 1 -0.4 
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-550 


3.3 ± 0.6 


2.8 +0.5 1 -0.4 


2.9 ± 0.4 
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-700 


0.9 ±0.1 


1.1 ±0.1 


1.0 ±0.1 
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-950 


0.28 ± 0.06 


0.23 4-0.05 / -0.04 


0.24 ± 0.04 


950 


-2700 


0.007 ± 0.003 


0.008 ± 0.004 


0.007 ± 0.003 



PT.fr 


[GeV] 


1/ff rf(7/JpT.ff [1/TeV] 








e±jets /i±jets 


<'±jets 





40 
170 


-40 

- 170 

- 1100 


14 ±2 14 ±2 14 ±2 
3.0 ±0.4 3.1 ±0.3 3.0 ±0.3 
0.050 ± 0.010 0.051 ±0.008 0.051 ± 0.008 









ila da jdyf, 






e±jets 


/i±jets 


l'±jets 




-2.5 --1 
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0.077 ± 0.009 


0.072 ± 0.008 


-1 --0.5 


0.32 ± 0.03 


0.35 ± 0.03 


0.34 ± 0.02 


-0.5 - 
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0.41 ± 0.02 


0.42 ± 0.02 
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0.43 ± 0.02 


0.42 ± 0.02 
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0.34 ± 0.03 


0.31 ± 0.02 


0.32 ± 0.02 


1-2.5 


0.080 ± 0.010 


0.083 ± 0.007 


0.080 ± 0.007 



Table 4 Con'elation coefficients between bins of the relative differential cross-section (top) l/cJ da /dnitj, (middle) \/a da /dpjjj and (bottom) 
l/a da jdyff in the combined I + jets channel. 
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-0.93 1.00 0.21 
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-0.40 0.24 
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Table 5 Percentage uncertainties on (top) I/ct d<y/dm,f, (middle) l/ff da/dpj ,f and (bottom) l/ff da /dy,j in the combined £'+jets channel. 
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